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We report on a search for the decay B° — > p°p° and other charmless modes with a tt + ty~tt + ty~ final 
state, including B° -> pV+vr", non-resonant B° -» 4^, B° -> p°/o(980), B° -> / (980)/ (980) 
and B° -> / (980)7r+7r . These results are obtained from a data sample containing 657 million BB 
pairs collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We set an 
upper limit on B(B° -> p°p°) of 1.0 x 10~ 6 at the 90% confidence level (C.L.). From our B° -> p°p° 
measurement and an isospin analysis, we determine the Cabibbo-Kobayashi-Maskawa phase <f>2 to 
be 91.7 ± 14.9 degrees. We find excesses in B° — > p°ir + ir~ and non-resonant B° — > 47T* with 1.3(7 
and 2.5cr significance, respectively. The corresponding branching fractions are less than 12.0 x 10 -6 
and 19.3 x 10~ 6 at the 90% C.L. In addition, we set 90% C.L. upper limits as follows: B(B° -> 
p°/o(980)) < 0.3 x 10~ 6 , B(B° -> f (980) f (980)) < 0.1 x 10~ 6 , and B(B° -> / (980)7r+7r-) < 
3.8 x 10" 6 . 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



In the Standard Model (SM), CP violation in the 
weak interaction can be described by an irreducible com- 
plex phase in the three-generation Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix Measure- 
ments of the differences between B and B meson decays 
provide an opportunity to determine the elements of the 
CKM matrix and thus test the SM. One can extract the 
CKM phase 2 = arg[-{V td V t * b )/(V ud V* b )] from the time- 
dependent CP asymmetry for the decay of a neutral B 
meson via ai-rn process into a CP eigenstate. How- 
ever, in addition to the b — > u process, there arc b — > d 
penguin transitions that shift the <p2 value by Sfa in the 
time-dependent CP violating parameter measurement. 
The shift Sfo can be determined from an isospin analysis 
[H of B — > 7T7T Q or B — > pp Q decays, or from a time- 
dependent Dalitz plot analysis of B — ► pir Q decays. 

For B — > pp decays, polarization measurements in 
B — > p + p~ 0] and i? ± — > / o ± p° @ show the dominance of 
longitudinal polarization, indicating that the final state 
in B — > p + p~ is very nearly a CP eigenstate. Measure- 
ments of the branching fraction, polarization and CP- 
violating parameters in B — > p°p° decays complete the 
isospin triangle. The tree contribution to B° — > p°p° is 
color-suppressed, so its branching fraction is expected to 
be much smaller than that for B — > or B 111 — > p^p - 

This also makes it especially sensitive to the penguin am- 
plitude, and using the B° 
an isospin analysis allows one to determine 
uncertainty from penguin contributions. 

Predictions for B° — > p /? using perturbative QCD 
(pQCD) Q or QCD factorization [2, 0] approaches sug- 
gest that the branching fraction B(B° — > p°p°) is at or 



(O ^ branching fraction in 
free of 



below 1 x 10~ 6 , and that its longitudinal polarization 
fraction /l is around 0.85. A non-zero branching frac- 
tion for B° —> p°P° has been reported by the BaBar 
collaboration [lOj, [ll]; they measured B(B° — > /O ^ ) = 
(0.92 ±0.32 ±0.14) x 10" 6 with a significance of 3.1 stan- 
dard deviations (a) , and a longitudinal polarization frac- 
tion, /l = 0.751^4 ± 0.05. They do not observe a non- 
resonant B° — > 47r ± or £>° — > p a ir + n~ contribution. The 
theoretical prediction for the non-resonant B a — > 47r ± 
branching fraction is around 1 x 10 -4 [l2j|- The most 
recent measurement of this decay was made by the DEL- 
PHI collaboration [3], which sets a 90% C.L. upper limit 
on the branching fraction of 2.3 x 10~ 4 . 

The data sample used in the analysis reported here 
contains 657 million BB pairs collected with the Belle 
detector at the KEKB asymmetric-energy e + e~ (3.5 and 
8 GeV) collider [l4|. o perating at the T(45) resonance. 
The Belle detector [15j,|l6| is a large-solid- angle magnetic 
spectrometer that consists of a silicon vertex detector, a 
50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel- like ar- 
rangement of timc-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to de- 
tect K° L mesons and to identify muons. Sig nal Monte 
Carlo (MC) is generated with EVTGEN |l7], in which 
final-state radiation is taken into account with the PHO- 
TOS package [l8j], and processed through a full detector 
simulation program based on GE ANT3 [3] . 

B° meson candidates are reconstructed from neutral 



3 



combinations of four charged pions. Charged track candi- 
dates are required to have a distance-of-closest-approach 
to the interaction point (IP) of less than 2 cm in the 
direction along the the positron beam (z-axis) and less 
than 0.1 cm in the transverse plane; they are also re- 
quired to have a transverse momentum px > 0.1 GeV/c 
in the laboratory frame. Charged pions are identified 
using particle identification (PID) information obtained 
from the CDC (dE/dx), the ACC and the TOF. We dis- 
tinguish charged kaons and pions using a likelihood ra- 
tio 7?.pid = £k /(£■!< + £tt), where C„(Ck) is the likeli- 
hood value for the pion (kaon) hypothesis. We require 
TipiD < 0.4 for the four charged pions. We require 
that charge tracks have a laboratory momentum in the 
range [0.5, 4.0] GeV/c, and a polar angle in the range 
[32.2, 127.2]°. For such tracks the pion identification effi- 
ciency is 90%, and the kaon misidentification probability 
is 12%. Charged particles that are positively identified 
as an electron or a muon are removed. 

To veto B —¥ D^n and B — > D s n backgrounds, we 
remove candidates that satisfy cither of the conditions 
IA/O^TttT) - m D(s) \ < 13 MeV/c 2 or [M^tt^) - 
m D o | < 13 McV/c 2 , where /i ± is either a pion or a kaon, 
and m-Dr,-) and m D o are the masses of the D( s ) and D° 
mesons, respectively. Furthermore, to reduce the B — > 
afn^ feeddown in the signal region, we require that the 
pion with the highest momentum have a momentum in 
the T(4S) center-of-mass (CM) frame within the range 
[1.30, 2.65] GeV/c. 

The signal event candidates are characterized by two 
kinematic variables: the beam-energy-constrained mass, 
A/ bc = yj -E^am — P"b i an d tne energy difference, AE = 
Eb — -Eboam, where -Ebcam is the run-dependent beam 
energy, and Pb and Eb are the momentum and en- 
ergy of the B candidate in the T(45) CM frame. In 
B° — > p°p° — > (7T + 7r~)(7r + 7r~) decays, or other charm- 
less modes with a 7r + 7r~7r + 7r~ final state, the invari- 
ant masses M(tt + it~) vs. M(ir + Tr~) are used to dis- 
tinguish different modes. There arc two possible com- 
binations for M(ir + ir~) vs. M(tt + tt~): (t^i^i)^} 71 ^) 
and (tti tt^ ) (vr^ ') , where the subscripts label the mo- 
mentum ordering, i.e. ir^(ir^) has a higher momentum 
than (ti2~). Here we consider both (itf Trf) (tt^tt^) 
and {^t 7T 2)( 7r 2 7r i ) combinations and select candidate 
events if either one of the combined masses lies within 
the signal window [0.55, 1.7] GeV/c 2 . This signal win- 
dow is chosen to accept p° — > 7r + 7r~, / (980) — > tt + tt~, 
and non-resonant modes, and to exclude K° — ► tt + tt~ 
and charm meson decays such as D° — » ir + ir~ . If both 
(7r^ 7r^)(7r^"7r^) and (^t 7T 2)( 7T 2 7T i) combinations of a 
candidate has a 7r + 7r~ pair with an invariant mass in 
the signal window, we select the (ir^ 7r^")(7r^7r^) com- 
bination. According to MC simulation, this criteria se- 
lects the correct combination for p°p° signal decays 98% 
of the time. For fitting, we symmetrize the M 2 (ir + ir~) 



vs. M 2 (tt + it~) Dalitz plot by plotting the tt^i ( 7T t 7T 2~) 
combination against the horizontal axis for events with 
an even (odd) event identification number, which is the 
location of the event in the data. 

The dominant background comes from continuum 
e + e~ — > qq (q = u, d, c or s) events. To distinguish 
signal from the jet-like continuum background, we use 
modified Fox- Wolfram moments (20| , which are combined 
into a Fisher discriminant. This discriminant is combined 
with PDFs for the cosine of the B flight direction in the 
CM frame and the distance in the z-axis between two B 
mesons to form a likelihood ratio 1Z = C S /{C S + £ q q). 
Here, C s (£ q q) is a likelihood function for signal (contin- 
uum) events that is obtained from the signal MC simula- 
tion (events in the sideband region M^ c < 5.26 GeV/c 2 ). 
We also use a flavor tagging qual ity variable r provided by 
the Belle tagging algorithm [21( that identifies the flavor 

of the accompanying B° meson in the T(45) — > B°B° 
decay. The variable r ranges from r = for no flavor 
discrimination to r — 1 for unambiguous flavor assign- 
ment, and it is used to divide the data sample into six 
r bins. Since the discrimination between signal and con- 
tinuum events depends on the r-bin, we impose different 
requirements on 1Z for each r-bin. We determine the 
1Z requirement such that it maximizes a figure-of-merit 
N s / s/N~ a + N q q, where N s (N q q) is the expected number 
of signal (continuum) events in the signal region. For 
22% of the events, there are multiple candidates; for these 
events we select the candidate with the smallest \ 2 value 
for the B° decay vertex reconstruction. This selects the 
correct combination 79.6% of the time. The detection 
efficiency for the signal is calculated by MC to be 9.16% 
(11.25%) for longitudinal (transverse) polarization. Since 
longitudinally polarized B° — > p° p° decays produce low 
momentum pions from one or both p s, their detection 
efficiency is lower than that for transversely polarized de- 
cays. 

Since there are large overlaps between B° — ► p°p° 
and other signal decay modes in the Mi(tt + it~) vs. 
M2(7r + 7r~) distribution, we distinguish these modes by 
fitting to a large M 1 (tt + tt~) vs. M 2 (Tr + Tr~) region. The 
signal yields are extracted by performing extended un- 
binned maximum likelihood (ML) fits. In the fits, we use 
four-dimensional (Mb C! AJ5, Mi, M^) information to dis- 
criminate among p°p°, p°Tr + ir~ , non-resonant 4-7T , p°fo, 
fofo and / 7r + 7r~ final states. We define the likelihood 
function 

\ j / i=1 \ j / 

where i is the event identifier, j indicates one of the 
event type categories for signals and backgrounds, rij 
denotes the yield of the j-th category, and Pj is the 
probability density function (PDF) for the j-th cate- 
gory. The PDFs are a product of two smoothed two- 
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dimensional functions: P] = Pj(M^ c , AE\ M\, Mj) = 
p(M^AE^xp(Ml,Aq). 

For the B decay components, the smoothed func- 
tions Psmoothod(M^ c , AE 1 ) and JWothedO^i, Mj) are ob- 
tained from MC simulations. For the M^ c and AE 
PDFs, possible differences between real data and the MC 
modeling are calibrated using a large control sample of 
B° — > L>~(A'+7r~7r~)7r+ decays. The signal mode PDF 
is divided into two parts: one is correctly reconstructed 
events and the other is "self-cross- feed" (SCF), in which 
at least one track from the signal decay is replaced by one 
from the accompanying B meson decay. We use different 
PDFs for the correctly reconstructed and SCF events, 
and fix the SCF fraction to that from the MC simulation 
in the nominal fit. 

For the continuum and charm B decay backgrounds, 
we use the product of a linear function for AE, an 
ARGUS function [22| for Mb c and a two-dimensional 
smoothed function for M1-M2. The parameters of the 
linear function and ARGUS function for the continuum 
events arc floated in the fit. Other parameters and the 
shape of the M1-M2 functions are obtained from MC 
simulations and fixed in the fit. 

For the charmless B decay backgrounds, we use three 
separate PDFs for B° — > afir T , B ± — » p^p and other 
charmless B decays; all of the PDFs are obtained from 
MC simulations. In the fit, we fix the branching fraction 
of B° -> af 71-T to the published value (33.2 ± 3.0 ± 3.8) x 
1(T 6 pj. If wc float the B° -> afn^ yield in the fit, 
the fit result is B{B a -> afir?) = (33.8+^) x 1(T 6 , 
which is consistent with the assumed value. Wc fix the 
yield of B ± - 

average branching fraction 
other charmless B decays. 

Tabic U and Fig. [T] show the fit results and projections 
of the data onto AE, M hc , Mi(7r+7r~) and M 2 (tt + tt~ ) for 
B° — > p a p° decay. The statistical significance is defined 
as y/— 2 ln(£o/£max), where Cq and £ max are the values 
of the likelihood function when the signal yield is fixed 
to zero and allowed to vary, respectively. The 90% C.L. 
(C.L.) upper limit for the yield N is calculated from the 
equation 



z p° to that expected based on the world 
and wc float the yield of 



Jo £{x)dx 
Jo°° £( x ) dx 



90%, 



(2) 



where x corresponds to the number of signal events. We 
include the systematic uncertainty into the upper limit 
(UL) by smearing the statistical likelihood function by 
a bifurcated Gaussian whose width is equal to the total 
systematic error. The significance including systematic 
uncertainties is calculated as before, except that we only 
include the additive systematic errors related to signal 
yield in the convoluted Gaussian width. 

The fractional systematic errors are summarized in Ta- 
ble |TTJ For the systematic uncertainties due to the fixed 



TABLE I: Fit results for the decay modes listed in the first 
column. The signal yields, reconstruction efficiencies (assum- 
ing the probability for the sub-decay mode /o(980) — > tt + tv~ 
is 100%), significance (S, in units of a), branching fractions 
(B, in units of 10 -6 ) and the upper limit at the 90% C.L. (UL, 
in units of 10 -6 ) are listed. For the yields and branching frac- 
tions, the first (second) error is statistical (systematic). 



Mode 



Yield 



Eff.(%) S B 



UL 



P°P° 

471"* 



9/ | r + 23. 6 + 10.1 
Z ^- J -22.1-16.2 



p u 7r + 7r- m.StesI ± 52 - 3 2 - 90 L3 ^-I'.l ± 2 -7 < 12.0 

161.2t^;^7.7 L9g 2 5 12 .4+4-7+2.1 < lg 3 

p°fo -ll-8t&S!t.8 9.81 - <0.3 

/o/o -7-7±tl ±3.0 10.17 - <0.1 

/ott+tt- 6.3±|X;? ± 18.0 2.98 - 0.3+\ ;| ± 0.9 < 3.8 




-0.05 0.05 

AE (GeV) 



5.26 5.28 

M bc (GeV/c 2 ) 




f.55 0.78 1.01 1.241.47 1.7 

M,(7tic) (GeV/c 2 ) 




f.55 0.78 1.01 1.241.47 1.7 

M 2 (jm) (GeV/c 2 ) 



FIG. 1: Projections of the four-dimensional fit onto (a) 
AE, (b) Mbc, (c) Mi(7r+7r~), and (d) M 2 (it + -k~), for candi- 
dates satisfying (except for the variable plotted) the criteria 
AE e [-0.05, 0.05] GeV, M bc G [5.27, 5.29] GeV/c 2 , and 
Afi, 2 (7r + 7r") 6 [0.626, 0.926] GeV/c 2 . The fit result is shown 
as the thick solid curve; the solid shaded region represents the 
B° — > p°p° signal component. The dotted, dot-dashed and 
dashed curves represent, respectively, the cumulative back- 
ground components from continuum processes, b — > c decays, 
and charmless B backgrounds. 



branching fractions, we vary the branching fractions of 
B° -> af 7T+ (33.2 ± 4A in units of 10" 6 ) [H and 
B ± -> p ± p° (18.2±3.0) H by their ±la errors. The fits 
are repeated and the differences between the results and 
the nominal fit values are taken as systematic errors. Sys- 
tematic uncertainties for the AE-Mhc PDFs used in the 
fit are estimated by performing the fits while varying the 
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signal peak positions and resolutions by ±lcr. System- 
atic uncertainties for the M1-M2 PDFs are estimated in 
a similar way. A systematic error for the longitudinal po- 
larization fraction of B° — > p°p° is obtained by changing 
the fraction from the nominal value ft, = 1 to the most 
extreme alternative value /l = 0. According to MC, the 
signal SCF fractions are 20.4% for B° -> p°p°, 14.2% 
for B° — > p°n + TT~, 11.1% for non-resonant B° — > 4vr , 
15.0% for B° -> p°f , 9.9% for B° -» f f and 13.4% for 
B° — > / 7r + 7r _ . We estimate a systematic uncertainty 
for the signal SCF by varying its fraction by ±50%. 

An MC study indicates that the fit biases are +2.4 
events for B° — > p°p°, +7.2 events for B Q — > p°7r + 7r+ 
+12.5 events for non-resonant B° — » 4tt ± , +3.6 events for 
£° -> p°/o, -0.8 event for B° -> / / and +5.1 events 
for i?° — > /o7r + 7r - . We find that fit biases occur due to 
the correlations between the two sets of variables (AE, 
Mbc) and (Mi, M2), which are not taken into account 
in our fit. We correct the fit yields for these biases. To 
take into account possible differences between the MC 
simulation and data, we take both the magnitude of the 
bias corrections and the uncertainty in the corrections as 
systematic errors. 

We study the possible interference between B° — > 
a^T: T , B° — > p°p°, B° — ► p°7r + 7r~ and non-resonant 
5° — ► 4tt ± using toy MC. We add a simple interference 
model to the toy MC generation, which is, for p° — > tt + it~ 
decay, modified from a relativistic Breit-Wigner function 
to 



1 



m 2 — mi + inioT 



Ae 



A 2 + 2A 



(to 2 — TOq) cos 5 — Tuiq sin 6 
(m 2 - ml) 2 + (r?7i ) 2 



(m 2 - m 2 ,) 2 + (rTO ) 2 



(3) 



where A and 6 arc the interfering amplitude and phase, 
and too and T are the p° mass and width, respectively. 
We assume that the interference term due to the ampli- 
tudes for B° — > af ir T , B° — > J o°7r + 7r~ and non-resonant 
B° — ► 47r ± decays are constant in the £?° — > p ^ sig- 
nal region. Since the magnitude of the interfering ampli- 
tude and relative phase are not known, we uniformly vary 
these parameters and perform a fit in each case to mea- 
sure the deviations from the incoherent case. We take 
the r.m.s. spread of the distribution of deviations as the 
systematic uncertainty due to interference. 

The systematic errors for the efficiency arise from the 
tracking efficiency, particle identification (PID) and TZ 
requirement. The systematic error on the track-finding 
efficiency is estimated to be 1.2% per track using par- 
tially reconstructed D* events. The systematic error due 
to the pion identification (PID) is 1.0% per track as es- 
timated using an inclusive D* control sample. The 1Z 



requirement systematic error is determined from the ef- 
ficiency difference between data and MC using a B° — > 
D~ (K + Tr~ir~)Tr + control sample. 



TABLE II: Summary of systematic errors (%) for the branch- 
ing fraction measurements, /l and /scf are the fractional 
uncertainties for longitudinal polarization and self-cross-feed. 



Source 




p°n + n 


471* 


p7o 


/o/o 


/o71 + 71~ 


Fitting PDF 


±10.2 


±29.8 


±12.2 


±18.6 


±31.2 


±270 


B(B° -» aiw) 


±21.6 


±33.5 


±2.7 


±17.8 


±1.3 


±39.7 


B{B ± -> pV) 


±0.0 


±0.7 


±0.2 


±0.0 


±0.0 


±1.6 


/l 


-53.7 












/scf 


±11.4 


±8.3 


±6.0 


±5.1 


±5.2 


±20.6 


Fit bias 


±16.3 


+6.4 
-5.7 


+7.8 
-3.3 


+30.5 
-14.4 


±20.8 


±82.5 


Interference 


+25.7 
-20.8 












Tracking 


±5.3 


±4.6 


±4.4 


±5.0 


±4.8 


±4.5 


PID 


±4.8 


±3.5 


±3.2 


±4.4 


±3.9 


±3.4 


1Z requirement 


±3.2 


±3.2 


±3.2 


±3.2 


±3.2 


±3.2 




±1.4 


±1.4 


±1.4 


±1.4 


±1.4 


±1.4 


Sum(%) 


+41.1 

-66.0 


±46.5 


+ 17.2 

-15.6 


+40.9 
-30.8 


±38.6 


±286 



To constrain 62 using B — » pp decays, we perform an 
isospin analysis 0, [25| using the measured branching frac- 
tions of longitudinally polarized B ± — ► p ± p°, B 
and B° — » p°p° decays as the lengths of the 
the isospin triangles. The B ± — > p ± p° and B 
branching fractions used, as well as the corresponding 
/l values, are world average values j24[; the B° — ► p°p° 
branching fraction is from this measurement, and we as- 
sume /l = I- The CP-violating parameters S~£~ and 



p+p 



p+p- 



are determined from the time evolution of the Ion 



gitudinally polarized B — > p + p~ decay [1, 0]. Fig. [2] 
plots the difference between one and the C.L. (1— C.L.) 
as a function of <p2', the central value and one sigma in- 
terval consistent with the SM is 2 = (91.7 ± 14.9)°. 
In summary, we measure the branching fraction of 



B° 



to be (0.4 ± 0.4 



+0.2 



1.3) x 



10~ b with l.Ocr sig- 



nificance; the 90% C.L. upper limit including systematic 
uncertainties is B(B° — ► p°p°) < 1.0 x 10" 6 . These val- 
ues correspond to longitudinal polarization (/l = 1); the 
upper limit is conservative as the efficiency for /l = 1 is 
smaller than that for /l = 0. If we take /l = 0.85, the 
average of the theoretical predictions 0, @] , the measured 
value becomes (0.3 ± 0.3) x 10~ 6 (statistical error only). 

On the other hand, we find excesses in B° — > /j°7r + 7r~ 
and non-resonant B° — » 47r ± decays with 1.3cr and 2.5a 
significance, respectively. We measure the branching 
fraction and 90% C.L. upper limit for B° — > /j°7r + 7r~ de- 
cay to be (5.9±|:|±2.7) x lO" 6 and B{B° -> p°Tr+Tr~) < 
12.0 x 10~ 6 . For the non-resonant B° — > 47r mode, we 
measure its branching fraction to be (12.4^4 g^ 2 'g) x 10 -6 
with a 90% C.L. upper limit of B(B° -> 4tt ± ) < 19.3 x 
10~ 6 . For these limits we assume the final state par- 
ticles are distributed uniformly in three- and four-body 
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100 120 140 160 180 
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FIG. 2: 1 — C.L. vs. <f>2(a) obtained from the isospin analysis 
of B — > pp decays. 



phase space. We find no significant signal for the decays 
B° -> p°fo, B° -> /o/o and B° -> /oTr+Tr"; the final 
results and upper limits are listed in Table |H 
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